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^ j , Abstract. In this paper we analyse five observations of the BL Lac object AO 0235+16 performed with the 
Chandra and XMM-Newton satellites during the years 2000-2005. In the February 2002 observation the source 

I is found in a bright state and presents a steep X-ray spectrum, while in all the other epochs it is faint and the 

^ . spectrum is hard. The soft X-ray spectrum appears to be strongly absorbed, likely by the intervening system at 

' z = 0.524, which also absorbs the optical-UV radiation. We find that models that consider spectral curvature are 

, superior to single power law ones in fitting the X-ray spectrum. In particular, we favour a double power law model, 

CO ' which agrees with the assumption of a superposition of two different components in the X-ray domain. Both in 

CO ' the Chandra and in one of the XMM-Newton observations, a tentative detection of the redshifted Fe Ka emission 

, line may suggest its origin from the inner part of an accretion disc. Thermal emission from this accretion disc 
might explain the UV-soft-X-ray bump that appears in the spectral energy distributions, when the X-ray spectra 

I are complemented with the optical-UV data from the Optical Monitor onboard XMM-Newton. More likely, the 

. bump can be interpreted in terms of an additional synchrotron component emitted from an inner region of the jet 

CL(' with respect to that where the lower-energy emission comes from. An inspection of the X-ray light curves reveals 

^ I that intraday variability occurs only when the source is in a bright state. 

H ■ 

' Key words, galaxies: active - galaxies: BL Lacertae objects: general - galaxies: BL Lacertae objects: individual: 

, AO 0235+16 - galaxies: jets - galaxies: quasars: general 



^ 1. Introduction 

rt , The BL Lac object AO 0235+16 (z = 0.94) is known to 
be an intriguing source, exhibiting a number of extreme 
properties. It has shown large amplitude variability on 
long time scales (months-years) and noticeable intraday 
variability (IDV) over a wide range of the electromagnetic 
spectrum, from the radio to the optical band (Takalo et al. 

Heidt & Wagner 1996: Noble & Miller IMF; Romero 
et al. .l997. Kraus et al. .l999. Romero et al. i2000i Raiteri 
et al. 120011 12005|l . These observations support the idea 
that the radiation from this object is strongly beamed to- 
wards us, with Dopplcr factors reaching up to ^ ~ 100 
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(Fujisawa et al. lT^ Frey et aL lMTUI Jorstad et al lMTTl 
Kraus et al OnUHjl . 

The interpretation of optical variability is made difh- 
cult by the presence of a complex environment (Yanny 
et al. iroS^ Burbidge et al. Nilsson et al. ITMSjl . 

A couple of objects at redshift z = 0.524 are very close 
to the line of sight: one is an AGN (named ELISA by 
Raiteri et al. I2UU5|I located 2 arcsec south of the source, 
affecting its optical-UV photometry. This object is prob- 
ably interacting with the other close object, a galaxy ly- 
ing 1.3 arcsec east of AO 0235+16, and this interven- 
ing system strongly absorbs the IR-to-soft-X-ray radia- 
tion of the blazar. Absorption lines due to this damped 
Lya (DLA) intervening system are seen in the optical and 
UV spectra of AO 0235+16. In particular, a fit to the 
Lya absorption line profile of the Hubble Space Telescope 
(HST) spectrum taken in February 1998 led to an esti- 
mate of the hydrogen column density at z = 0.524, which 
is A'^Hi = 4.5 ± 0.4 X 10^^ cm^^ according to Turnshek 



C. M. Raiteri et al.; X-ray emission from AO 0235+16: XMM-Newton and Chandra 



2 

et al. (12003), and iVni = 5 ± 1 x 10^^ cm"^ according to 
Junkkarinen et al. 1)2004(1 . 

Moreover, objects in the same intervening system 
could act as microlenses for the source emission, so that 
an interpretation of at least some episodes of variability in 
terms of extrinsic processes cannot be ruled out (Stickel 
et al. [TOSSI Takalo et al. [TUM Webb et al. HUM but see 
also the critical discussion by Kavser [T988|l . 

AO 0235+16 has recently been the target of a huge 
multiwavelength observing effort led by the Whole Earth 
Blazar Telescope^ (WEBT; e.g. Villata et al. HUM E0021 
[2004a, 2004bl) in the observing seasons 2003-2004 and 
2004-2005. Three XMM-Newton pointings were awarded 
to get information on the high-energy emission of the 
source. The results of the radio to optical monitoring dur- 
ing the first season of the campaign, together with an anal- 
ysis of the XMM-Newton pointing of January 2004, have 
already been published by Raiteri et al. ( 2005). Here we 
present the results of the XMM-Newton observations per- 
formed in August 2004 and January 2005, during the sec- 
ond observing season of the WEBT campaign. The other 
multiwavelength data collected by the WEBT in this pe- 
riod will be published elsewhere (Raiteri et al. I2006|l . We 
also analyse an archival XMM-Newton pointing performed 
in February 2002 and re-analyse the Chandra data taken 
in August 2000 (Turnshek et al. 12003(1 and those from the 
XMM-Newton pointing of January 2004. The aim is to 
better understand both the X-ray properties of this enig- 
matic source through a homogeneous analysis of the most 
recent data and their relationship with the lower-energy 
emission. In particular, we investigate the existence of a 
spectral curvature in the X-ray energy range as well as 
the presence of a bump in the UV-soft-X-ray region of 
the SED. 

This paper is organized as follows: in Section 2 we re- 
view past observations of AO 0235+16 by X-ray satellites. 
A detailed analysis of the X-ray data acquired by Chandra 
in August 2000 and by XMM-Newton in 2002-2005 is per- 
formed in Section 3, where the possibility that the fluores- 
cent Fe Ka line was detected in some of the above spectra 
is also discussed. In Section 4 the X-ray spectra are com- 
plemented by the optical and UV data from the Optical 
Monitor onboard XMM-Newton and the resulting UV- 
soft-X-ray bump is discussed in Section 5. The X-ray light 
curves corresponding to the above pointings are presented 
in Section 6. Conclusions are drawn in Section 7. 



2. Previous X-ray observations by otiier satellites 

In this section we briefly review the results derived from 
past X-ray observations of AO 0235+16. They are summa- 
rized in Tabled in terms of the hydrogen column density 
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A^H adopted to model the X-ray absorption, the photon 
index F^, and the unabsorbed 1 keV flux density. 

Worral & Wilkes ((1990(1 analysed the data taken by the 
Einstein Observatory imaging proportional counter (IPC) 
in 1979-1981. They fit the data with a power law with free 
Galactic absorption and noticed that an A'h value larger 
than the typical one for Galactic absorption was required 
to obtain a good fit. 

EXOSAT observed the source in 1983-1986. Data 
taken during the pointing of August 1, 1984 were anal- 
ysed by Ghosh & Soundararajaperumal 1(1995(1 in terms 
of a simple power law plus absorption; the highly uncer- 
tain value of A^H was compatible with the Galactic one. 

Madejski et al. (,1996 ) presented data taken by ROSAT 
in summer 1993 and data obtained by ASCA in February 
1994. They first fit the X-ray spectrum with a power 
law plus free absorption model, and this led to an extra- 
absorption. Then, a model with fixed Galactic absorp- 
tion plus free absorption at z = 0.524 was adopted. The 
ROSAT data showed AO 0235+16 in a bright and ac- 
tive state, revealing flux variations up to ~ 80% on a 3 
day timescale. On the contrary, no short-term variability 
was detected by ASCA, which observed the source in a 
much fainter state. The comparison between the ROSAT 
and ASCA results revealed spectral variability on the long 
term, which was interpreted by Madejski et al. ((1996(1 as 
a dominance of the steep tail of the synchrotron compo- 
nent when the source is bright, and a prevalence of the 
fainter and harder Compton component when it is faint. 
The authors ruled out the possibility that a significant Fe 
K line was present at either z = 0.524 or z = 0.94, with 
90% confidence limits on the line equivalent width of ^ 60 
eV and ~ 100 eV, respectively. 

Because of the flux variability detected in the ROSAT 
data, Comastri et al. ^2S7j analysed the ROSAT observa- 
tions of summer 1993 by separating those with low count 
rates from those with high count rates. Indeed, they re- 
ported a factor 3 difference in the flux, but a more or 
less constant spectral index. 

The source was observed by RXTE in November 1997, 
during the rising phase of a big radio-optical outburst 
(Webb et al.EOOOl Raiteri et al. OnnT)l . Webb et al. (PHT))! 
analysed the X-ray data with a power law model and re- 
ported a very low X-ray flux (see Table^l . It was suggested 
that the X-ray response to the low-energy fiux increase 
may be delayed or insignificant. 

ASCA pointed again at the source on February 11, 
1998; the results were published by Junkkarinen et al. 
(120041 . The spectrum was harder with respect to the pre- 
vious ASCA pointing, 4 years earlier, and the flux in the 
2-10 keV energy range was 3.4 times higher. 

Data acquired by BeppoSAX on January 28, 1999 were 
reported by Padovani et al. (2004). They found extra- 
absorption and both spectral index and flux consistent 

^ The photon index F is related to the energy index a 
through r = Of + 1, recaUing that the specific flux at frequency 
u is (X u~°' 
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Table 1. Results from previous X-ray observations of AO 0235-1-16 by other satellites. When two numbers appear 
in the A'h column, the first refers to the Galactic absorption, the second to absorption by the z — 0.524 intervening 
system. Fx j^ev is the unabsorbed 1 keV flux density 



Satellite 


Obs6rv£Ltioii 


Energy range 




r 


-^1 kcV 








fkeV") 
V ; 






ill Tvl 




Einstein 


1979-1981 


0.1-3.5 


> 0.76 


3.25+t« 


1.51 


(1) 


EXOSAT 


Aug. 1, 1984 


0.1-10 


8t? 


1 yr+O.SO 


i-44lt:^^ 


(2) 


ROSAT 


Jul. 21 - Aug. 15, 1993 


0.1-2.1 


0.76, 3.8+i i 


2.67±0.32 


1.24 


(3) 


ROSAT 


Jul. 26, 1993 


0.1-2.0 


2.46t;«r 


Z.03_Q 7g 


1.15 


(4) 


ROSAT 


Aug. 15, 1993 


0.1-2.0 


2.6it?:l 


2 =4+0.88 


0.41 


(4) 


ASCA 


Feb. 4-19, 1994 


0.5-10 


0.76, 4.6±1.3 


1.96±0.09 


0.3 


(3) 


RXTE 


Nov. 3-6, 1997 


2-60 




1.41±0.20 


(0.14±0.02)'' 


(5) 


ASCA 


Feb. 11-12, 1998 


0.5-10 


0.76, 4.6±1.3 


1.72±0.07 


0.278 


(6) 


BeppoSAX 


Jan. 28, 1999 


0.1-10 


2.8 






(7) 



" (1) Worral fc Wilkes (19901: (2'1 Ghosh fc Soundararaiaperumal (19951: (3) Madeiski et al. II1996II : 

(4) Comastri et al. (il997il : (5) Webb et al. (i2000J : (6) Junkkarinen et al. (.2004J : (7) Padovani et al. (IMTH 
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with the values derived by Madejski et al. (|1996|l for the 
ASCA data. The problem of a possible contribution to the 
X-ray flux by the southern AGN was also discussed, lead- 
ing to the conclusion that it should be roughly a factor 10 
smaller than that observed for AO 0235+16. 

We notice that Tabled can help to trace both the flux 
and spectral X-ray variability of the source, but it is not 
exhaustive, since in most of the above mentioned papers 
more than one model was applied and discussed. Yet, it 
is not possible to find a common model for all cases, and 
the comparison among results from different observations 
suffers from inhomogeneous data analysis. It is to avoid 
such a difficulty that in this paper we have analysed or re- 
analysed the last five X-ray observations using the same 
data treatment'^. 



3. XMM-Newton and Chandra X-ray spectra 

In this section we present the X-ray spectra of AO 
0235+16 taken by XMM-Newton in August 2, 2004 and 
January 28, 2005, during the second observing season of 
the WEBT campaign on this source. The results will be 
compared with those obtained by re-analysing the XMM- 
Newton observation performed in January 18-19, 2004, 
during the first season of the WEBT campaign (Raiteri 
et al. I2005|l . and with those derived from the analysis 
of archival XMM-Newton spectra taken in February 10, 
2002. Moreover, we also re-analyse the data of a Chandra 
pointing on August 20-21, 2000, which have been already 
presented by Turnshek et al. 112003). Table [21 reports the 
date, the UT of start and end as well as the duration of 
the exposure for each pointing and for each instrument. 



^ A uniform re-analysis of the archival ASCA and BeppoSAX 
data is given in Kadler (120051 



Table 2. Log of the Chandra and XMM-Newton obser- 
vations 

Detector Start End Duration 

(UT) (UT) (s) 

Chandra: August 2000 

PI: D. A. Turnshek 

ACIS-S Aug. 20, 23:39:52 Aug. 21, 08:43:05 30625 ~ 
XMM-Newton: February 2002 
PI: M. Watson 

MOSl Feb. 10, 04:27:48 Feb. 10, 09:46:33 19125 
MOS2 Feb. 10, 04:27:47 Feb. 10, 09:46:33 19126 

pn Feb. 10, 05:01:03 Feb. 10, 09:41:50 16847 

XMM-Newton: January 2004 
PI: C. M. Raiteri 
MOSl Jan. 18, 19:07:18 Jan. 19, 03:25:01 29863 
MOS2 Jan. 18, 19:07:18 Jan. 19, 03:25:06 29868 
pn Jan. 18, 19:12:28 Jan. 19, 03:26:59 29671 

XMM-Newton: August 2004 
PI: C. M. Raiteri 
MOSl Aug. 02, 00:42:16 Aug. 02, 03:56:35 11659 
MOS2 Aug. 02, 00:42:14 Aug. 02, 03:56:40 11666 
pn Aug. 02, 00:47:24 Aug. 02, 03:58:35 11471 

XMM-Newton: January 2005 
PI C. M. Raiteri 
MOSl Jan. 28, 14:43:36 Jan. 28, 19:17:57 16461 
MOS2 Jan. 28, 14:43:40 Jan. 28, 19:18:03 16463 
pn Jan. 28, 14:48:46 Jan. 28, 19:19:56 16270 



3.1. Data reduction 

Chandra observed AO 0235+16 on August 20-21, 2000 
with ACIS-S. The image, which has a space resolution of 
0.492 arcsec/px, docs not show any emission 2 arcsec south 
of the source, thus ruling out an appreciable contribution 
from the southern AGN, named ELISA. 

Data were reduced with the CIAO package version 
3.0.2. We extracted the source spectrum from a circular 
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region of 3 arcsec size around the source and the back- 
ground spectrum in a surrounding annulus of radii 5 and 
10 arcsec. Only the energy channels from 0.4 to 8.0 keV 
were considered. 

The European Photon Imaging Camera (EPIC) on- 
board XMM-Newton has three detectors: MOSl, M0S2, 
and pn; they all point at the same object at the same time. 
All XMM-Newton observations analysed in this paper but 
the one of February 10, 2002, were performed in small- 
frame, medium-filter configuration in order to avoid pos- 
sible problems of pile-up^, out-of-time events, and degra- 
dation. 

The data were reduced with the Science Analysis 
System (SAS) software version 6.0 ^. Source spectra were 
extracted from a circular region centered on the source; 
the aperture radius was 30-35 arcsec in all cases but for 
the February 2002 observations, where it was 100 arcsec. 
Background spectra were derived from source-free regions 
on the same CCD of the source. We considered only chan- 
nels with energy greater than 0.2-0.3 keV and lower than 
10-12 keV, depending on the instrument (MOS having 
stricter limits than pn) and data quality. 

In the February 10, 2002 observation the source was 
detected in a high X-ray state. As the observation was 
performed using a full-frame, thin-filter configuration, we 
also checked the possible presence of pile-up with the SAS 
task epatplot. We found a small excess of double events in 
the ~ 1.5-5 keV energy range in the pn observation, which 
was removed by extracting the source spectrum in an an- 
nulus with an inner radius of 5 arcsec. However, the pn 
spectrum appears somewhat distorted below ~ 0.5 keV, 
hence for the spectral analysis of this observation we used 
the data from MOSl and M0S2 only. 

After this preliminary data reduction with the specific 
packages, each source spectrum was associated with its re- 
sponse matrix file, ancillary response file, and background 
spectrum, as well as binned in order to have a minimum 
of 25 counts in each bin. This was accomplished with the 
task grppha of the FTOOL package. These "grouped" spec- 
tra were then analysed with the Xspec package, version 
11.3.0. For each epoch, only the reliable energy channels 
were selected, discarding measurements affected by large 
uncertainties. For all the XMM-Newton observations, data 
from the different EPIC detectors are fitted simultane- 
ously, i.e. the same model is apphed to the MOSl, M0S2, 
and pn spectra to increase the statistics. 

The resulting spectra are shown in Figs. ^El fitted by 
a double power law kiE~^^ + k2E~^^ with Galactic ab- 
sorption A'h = 0.679 X 10^^ cm~^ plus extra-absorption by 

* Photon pile-up occurs when more than one X-ray photon 
arrives in one camera pixel or in adjacent pixels before it is 
read out; it can affect both the source point spread function 
and the spectral response of the instrument. 

^ When extracting the spectrum of both the source and the 
background, the strings "FLAG==0" and "PATTERN<=4" were 
included in the selection expression for the evselect task of 
the SAS software for all the three EPIC detectors to provide 
conservative event selection criteria. 



Chandra: August 20, 2000 

ACI5 




Fig. 1. The X-ray spectrum of AO 0235+16 detected by 
Chandra on August 20-21, 2000 fitted by a double power 
law model with fixed Galactic and z = 0.524 DLA system 
absorptions; the lower panel shows the ratio of the data 
to the model fit 
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Fig. 2. The X-ray spectrum of AO 0235+16 detected by 
XMM-Newton on February 10, 2002 fitted by a double 
power law model with fixed Galactic and z = 0.524 DLA 
system absorptions; the lower panel shows the ratio of the 
data to the model fit 

the intervening system at z = 0.524 — 5.0 x 10^^ cm~^ 
(see next section). The ratio of the data to the model fit 
is shown in the bottom panels. 

In Figs.lSElthe MOSl and M0S2 spectra are plotted 
with black squares and red triangles, respectively, while 
the pn one is shown with green diamonds. 

3.2. Results 

We start by considering 

— Model 1: a single power law kE~^ 

with fixed Galactic photo-electric absorption, and extra- 
absorption attributed to the intervening DLA system at 
z = 0.524. We adopted the value iVn = 0.679 x lO^^ cm^^ 
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XMM-Newton: Jonuory 18-19. 2004 
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Fig. 3. The X-ray spectrum of AO 0235-M6 detected by 
XMM-Newton on January 18-19, 2004 fitted by a double 
power law model with fixed Galactic and z = 0.524 DLA 
system absorptions; the lower panel shows the ratio of the 
data to the model fit 



XMM-Newton: August 2, 2004 
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Fig. 4. The X-ray spectrum of AO 0235-1-16 detected by 
XMM-Newton on August 2, 2004 fitted by a double power 
law model with fixed Galactic and z = 0.524 DLA system 
absorptions; the lower panel shows the ratio of the data 
to the model fit 

for the Galactic hydrogen column, from the Leiden- 
Argentinc-Bonn Survey^ (Kalberla et al. i2()05|l . which is 
the most sensitive Milky Way Hi survey to date, with the 
most extensive coverage both spatially and kinematically. 
Both the Galactic and the z = 0.524 absorber are first as- 
sumed to have solar abundance ratios, as was done in the 
past data reductions which led to the results summarized 
in Tabled 

Table 131 shows the results for the Chandra and XMM- 
Newton pointings analysed in the present work (first line) . 
The fit is good for all cases, but the best-fit A'h value of 
the z = 0.524 system is nearly half the A^hi value derived 
by fitting the Lya absorption line in the HST spectrum of 

^ http : //www. as tro ■ uiii-boim.de/~webrai/gerineLii/tools I 
_labsurvey . php 



M0S1+U052+pn 




Fig. 5. The X-ray spectrum of AO 0235+16 detected by 
XMM-Newton on January 28, 2005 fitted by a double 
power law model with fixed Galactic and z — 0.524 DLA 
system absorptions; the lower panel shows the ratio of the 
data to the model fit 



February 1998, which is 4.5 ±0.4 x 10^^ cm~^ according to 
Turnshek et al. H2003|l . and 5± 1 x 10^^ cm~^ according to 
Junkkarinen et al. H20()4|l . Notice that the A^hi value de- 
rived from the Lya absorption line accounts for neutral hy- 
drogen absorption only, and it is thus a lower limit to the 
A^H value derived with Xspec, which instead takes into ac- 
count all the possible absorptions of all elements in molec- 
ular, neutral, or ionized form. Moreover, the Galactic gas 
is metal-poor with respect to solar, and a metal-poor 
chemical composition appears more suitable also to de- 
scribe the z = 0.524 absorber. We thus re-analyse our 
datasets by using the Wilms et al. POOO) model for both 
the interstellar medium and the z = 0.524 absorber. This 
model includes, besides revised abundances for the inter- 
stellar medium, updated photoionization cross-sections, as 
well as a treatment of interstellar grains and H2 molecules. 
For each observation in Tabled the results of the corre- 
sponding analysis are shown in the second row. As one can 
see, there are little differences in the spectral indices and 
fluxes with respect to the previous case, but the value of 
A^H increases, as expected, since now we have less absorb- 
ing metals. However, it is still well below the value of A^hi 
derived by Turnshek et al. (20(3|l and Junkkarinen et al. 
(2004)7. 

We then fix the parameter N^i^"^^ = 5.0 x 10^^ cm~-^. 
As a consequence, the single power law model gives a worse 
fit (see third rows in Table El), especially for the XMM- 
Newton spectra of February 2002. Thus, we investigate 
the reliability of other models implying some spectral cur- 
vature. We consider 

- Model 2: a logarithmic parabola fc£:(-r+/3iogi5) 
(Landau et al. 119861 Giommi et al. IMl^ Perlman et 
al. .200511 

7 When we try more complicated models introducing some 
spectral curvature, the A^h values increase a bit more, but the 
uncertainties of the parameters are large. 
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— Model 3: a broken power law kE for E < inbreak, 

kEll'l'E-^' for E > inbreak 

— Model 4: a curved model fci^-^i (1 + i;/£'brcak)^' 
(Fossati et al. I^ITOl Tavecchio et al. IM1T|) 

— Model 5: a double power law kiE~^^ + k2E~^^ 

all of them with the same two absorption components of 
Model 1, with fixed A'h values and with the Wilms et al. 
H2Ul)U|l model for the photo-electric absorption. The results 
of these fits are shown in Tables 01[7| The last column 
reports the F-test probability of the model with respect 
to the Model 1 case with N^^"^^ = 5.0 x lO^^ cm-^ (third 
rows in Table OJ. The lower the probability, the better the 
new model is with respect to the previous one. As one can 
see from the statistics, all these models fit the data better 
than a single power law. However, there is not a unique, 
best-fit model for all the available data. 

We notice that this spectral curvature is expected 
from the superposition of the high-energy tail of the 
synchrotron (or thermal, see Sect. 6) component on the 
(higher-energy) inverse-Compton branch, both compo- 
nents possibly implying power-law spectral shapes in this 
limited energy range. Moreover, in the specific case of AO 
0235-1-16 the curvature should be more dramatic, since 
previous works have highlighted the existence of a UV- 
soft-X-ray bump in the SED, indicating that the soft X- 
ray spectrum must turn up very abruptly towards the 
UV spectral region (Junkkarinen et al. 12(7041 Raiteri et 
al. I2005|l . From this point of view. Model 5, being the 
sum of two power laws, appears to be the most suitable 
model to describe the source spectrum. 

The results of our study are also presented in Figs. 
Hil- llOl where the model spectra are shown for each obser- 
vation. Different model fits are plotted in different colours: 
Model 1 (single power law) is represented in red. Model 
2 (logarithmic parabola) in cyan, Model 3 (broken power 
law) in green. Model 4 (curved) in pink, and Model 5 (dou- 
ble power law) in blue. Unfolded and unabsorbed data ob- 
tained by applying the double power law model are also 
plotted, after some binning; we have to warn that had we 
chosen another model to unfold the data, the points would 
have shown accordingly. 

3.3. Comparison with other works 

By comparing the results of Table |31 with those of Tabled 
we notice that the two ASCA pointings of February 1994 
and 1998 detected the source at fiux levels comparable to 
those observed by XMM-Newton in January and August 
2004 (Fikcv ^ 0.3/iJy). The ASCA data were analyzed 
by Madejski et al. (T^) and Junkkarinen et al. H2004|l 
with prescriptions for absorption similar to ours, hence 
the results can easily be compared, highlighting how the 
source can exhibit noticeable spectral variations (from F ^ 
1.6 to F ~ 2.0) at the same brightness level. 

A previous analysis of the August 2000 Chandra ob- 
servation was published by Turnshek et al. (^^Hl- The 
interest of these authors was to study the metallicity of 
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Fig. 6. Fits to the X-ray spectrum of AO 0235-1-16 during 
the Chandra observation of August 20-21, 2000 accord- 
ing to Models 1-5; the blue dots represent the unfolded 
and unabsorbed data obtained by applying Model 5, after 
some binning 



XMM-Newton: February 10, 300S 




Fig. 7. Fits to the X-ray spectrum of AO 0235-1-16 dur- 
ing the XMM-Newton observation of February 10, 2002 
according to Models 1-5; the blue dots represent the un- 
folded and unabsorbed data obtained by applying Model 
5, after some binning 



the intervening DLA system dX z — 0.524. They fit the 
Chandra spectrum with a power law in the 2-8 keV en- 
ergy range, where absorption is expected not to affect the 
continuum flux, in order to flx both the slope and nor- 
malization of the X-ray spectrum. A best-fit photon in- 
dex F = 1.61 ± 0.10 was found. They adopted a Galactic 
Hi column density of 0.87 x 10^^ cm~^ after Hartmann & 
Burton ifTMTjl and a value of A'^hi = 4.5 ± 0.4 x lO^^ cm'^ 
for the z — 0.524 absorber, which was derived from a fit of 
the Lya absorption line in the February 1998 HST spec- 
trum. Then they tried different prescriptions for the metal 
content of both the Galaxy and DLA system, finding a 
range of possible metallicity for this latter of 10-40% of 
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Table 3. Results of the X-ray data analysis with Model 1: a single power law with Galactic absorption of 
0.679 X 10^^ enr^^ and extra-absorption at 2: = 0.524; for each observation the first row is obtained by adopting 
solar abundances, while the second row refers to the metal-poor model by Wilms et al. H2000|l . the third row fixing 
the extra-absorption at 5.0 x 10^^ cm~^. Fi^ev is the unabsorbed 1 keV fiux density; F is the observed 2-10 keV flux. 
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Table 4. Results of the X-ray data analysis with Model 2: a logarithmic parabola with absorption according to the 
metal-poor model by Wilms et al. (J2000j). The Galactic hydrogen column is 0.679 x 10^^ cm^^, while that of the 
z — 0.524 absorber is 5.0 x 10^^ cm^^. i^ikcV is the unabsorbed 1 keV flux density; F is the observed 2-10 keV flux. 
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Table 5. Results of the X-ray data analysis with Model 3: a broken power law with absorption according to the metal- 
poor model by Wilms et al. (2000). The Galactic hydrogen column is 0.679 x 10^^ cm^^, while that of the z = 0.524 
absorber is 5.0 x lO'^^ cm~^. FikcV is the unabsorbed 1 keV flux density; F is the observed 2-10 keV flux. 
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the solar one. The results of Turnshek et al. I)2()03|l are in 
fair agreement with those we found by applying a similar 
model (single power law with iV^-^^"* = 5.0 x lO^^cm"^, 
see Table EJ. 

Foschini et al. H2006|l analysed the data of the XMM- 
Newton observation of February 2002. When fitting a sin- 
gle power law with Galactic absorption in the 1-10 keV 
energy range, they obtained F = 2.28 ±0.02, roughly con- 
sistent with the results reported in Table |2| 



3.4. Iron lines 

The detection of Ka iron fiuorescent emission lines in 
AGN X-ray spectra is a very important and unique tool 
to study the immediate vicinity of the supermassive black 
hole which most likely is the central engine of the AGN 
itself (see e.g. the review papers by Fabian et al. i2000i and 
Reynolds & Nowak l^UnSjl . Up to now, Ka iron fluorescent 
lines have been detected in several Seyfert galaxies (see 
the above reviews and references therein) and in some ra- 
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Table 6. Results of the X-ray data analysis with Model 4: a curved model with absorption according to the metal-poor 
model by Wilms et al. 1)2000(1 . The Galactic hydrogen column is 0.679 x 10^^ cm^^, while that of the z — 0.524 absorber 
is 5.0 X 10^^ cm~^. i^ikcV is the unabsorbed 1 keV flux density; F is the observed 2-10 keV flux. 
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Table 7. Results of the X-ray data analysis with Model 5: a double power law with absorption according to the metal- 
poor model by Wilms et al. (2000). The Galactic hydrogen column is 0.679 x 10^^ cm~^, while that of the z = 0.524 
absorber is 5.0 x 10^^ cm^^. FikcV is the unabsorbed 1 keV flux density; F is the observed 2-10 keV flux. 
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XMM-Newton: January 18-19, 2004 
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Fig. 8. Fits to the X-ray spectrum of AO 0235-1-16 dur- 
ing the XMM-Newton observation of January 18-19, 2004 
according to Models 1-5; the blue dots represent the un- 
folded and unabsorbed data obtained by applying Model 
5, after some binning 
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XMM-Newton: August 2, 2004 
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Fig. 9. Fits to the X-ray spectrum of AO 0235-1-16 during 
the XMM-Newton observation of August 2, 2004 accord- 
ing to Models 1-5; the blue dots represent the unfolded 
and unabsorbed data obtained by applying Model 5, after 
some binning 



dio galaxies (e.g. Eracleous et al. I2000p . As far as we are 
aware, among blazars, only 3C 273 has shown this line in 
its X-ray spectrum (e.g. Yaqoob & Serlemitsos 2000), and 
it has been never seen in BL Lac objects, except for the 
possible case of S5 0716-K71 (Kadler et al. lMT^ . If BL Lac 
objects have an accretion disc surrounded by a hot corona 
and if this disc is cold enough to allow the emission of 
Fe Ka fluorescent line, then this line could be hidden by 
the beamed X-ray continuum of the jet (of synchrotron or 
inverse-Compton nature), just like the softer (IR-to-UV) 



emission lines are swamped by the corresponding jet con- 
tinuum. Alternatively, the disc in BL Lacs, if any, could 
be hot, or optically thin, or highly ionized, thus not pro- 
ducing the Ka feature. 

AO 0235-f 16 has been classified as a BL Lac object. 
However, it has shown both narrow and broad emission 
lines in its optical spectrum (Cohen et al. 119871 Nilsson 
et al. 1996; Raiteri et al. l2006|l . in particular a prominent 
Mgii broad line. These features are obviously better seen 
when the optical continuum from the jet is lower. The 
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XMM-Newton: January 28, 2005 
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Fig. 10. Fits to the X-ray spectrum of AO 0235+16 dur- 
ing the XMM-Newton observation of January 28, 2005 
according to Models 1-5; the blue dots represent the un- 
folded and unabsorbed data obtained by applying Model 
5, after some binning 



four X-ray spectra taken with Chandra in 2000 and with 
XMM-Newton in 2004-2005 show the source in a very 
faint X-ray state. In both the Chandra and the August 
2004 XMM-Newton spectra (see Figs. Q] and a possible 
emission feature is present around 2.8 keV (and a more 
marginal one around 3.1 keV). By adding a Gaussian line 
of 0.1 keV width in the energy interval 2.7-3.2 keV to 
Model 1 (case with N^^^^ = 5.0 x 10^^ cm"^)^ .^he best fit 
is reached when the line is at 2.8 keV and implies some im- 
provement of x^- — 2.63 for the Chandra observation, 
and 2.25 in the case of the August 2004 XMM-Newton 
one. The corresponding equivalent widths are ~ 61 eV 
and ^ 35 eV, respectively, but with large uncertainties. 
Consequently, we can only say that a line may be present 
in these two observations, with EW less than ~ 120 eV 
and ~ 75 eV, respectively. This possible emission line in 
the source frame would have an energy of about 5.4 keV, 
a bit far from the rest-frame energy of 6.40 keV of the cold 
Fe Ka fluorescent line. However, it is well known that, if 
the iron lines are produced in the inner accretion disc, they 
are affected from both Doppler and gravitational relativis- 
tic effects. Indeed, although they are intrinsically narrow, 
they can be broadened and skewed by the large orbital 
velocities of the accreting material, and redshifted by the 
strong gravitation field of the black hole. In fact, iron lines 
have been detected peaking at energies from 6-7 keV down 
to well below 6 keV in the source frame, also with a va- 
riety of shapes and widths (see the review by Fabian et 
al. I2UU0I and references therein). Even in the same source, 
the detected line can show dramatic changes in flux, pro- 
file, and redshift with time, as in the emblematic case of 
the Seyfert 1 galaxy MCG —6-30-15 (e.g. Iwasawa et al. 

mm . 



4. Ultraviolet to X-ray spectral energy distribution 

Besides the X-ray detectors, the XMM-Newton satellite 
also carries an Optical Monitor (OM), which allows to 
obtain photometric data in the optical U, B, V and ultra- 
violet UVVFl, UVM2, \JYW2 bands simultaneously with 
the X-ray observations. Optical and UV grisms are also 
available. 

The optical and ultraviolet emission from AO 0235+16 
is contaminated by the emission of ELISA, the AGN which 
is located only about 2 arcsec southward of the source. The 
ELISA contribution becomes more important as long as 
AO 0235+16 gets fainter, and as long as higher frequen- 
cies are considered (see discussion in Raiteri et al. I2005|l . 
The XMM-Newton source magnitudes listed in the OM 
pipeline products are obtained by performing aperture 
photometry with a 12 pixel (about 5.7 arcsec) aperture 
radius, so that they include the contribution of ELISA. 
In order to measure the magnitude of AO 0235+16 avoid- 
ing the southern AGN contamination, we analysed the 
OM images by following the same procedure adopted in 
Raiteri et al. (2005), i.e. by using the IRAF routine phot 
with an aperture radius of 2 pixels. We performed differen- 
tial photometry with respect to Stars 1, 2, 3 (Smith et al. 
CD85), and 6 (Fiorucci et al. ll998|l . using for these stars the 
standard magnitudes listed in the pipeline products. The 
source magnitude was derived by averaging the four values 
obtained with respect to the reference stars. Table |H| re- 
ports the inferred optical-UV magnitudes of AO 0235+16 
during the XMM-Newton pointings. 

The ultraviolet to X-ray spectral energy distributions 
are shown in Fig.^] The OM magnitudes of Table|S|have 
been dereddened according to Junkkarinen et al. H2004|l . 
using the extinction values listed in Table 5, Col. 5 of 
Raiteri et al. H2()()5(l . Then the OM fluxes have been ob- 
tained with respect to Vega. Notice that in Fig. ^2 we 
have lowered the UVA'/2 flux by 10% to take into ac- 
count the fact that at this frequency the Junkkarinen et 
al. H2()()4|l prescriptions overestimate the extinction. The 
X-ray spectra are shown according to Model 5 (double 
power law) , the spectrum thickness representing the stan- 
dard deviation of the data from the model fit, which is 
(7 ~ 0.049 for the Chandra pointing of August 2000, and 
cr = 0.024, 0.033, 0.046, and 0.065 for the XMM-Newton 
pointings of February 2002, January 2004, August 2004, 
and January 2005, respectively. 

One can see how the optical-UV state is higher when 
the X-ray state is higher. Moreover, the January and 
August 2004 data show how the steep optical spectrum has 
a minimum in the U-VYWl frequency range and then the 
spectrum becomes hard in the UV band, which is also con- 
firmed by the January 2005 data. In other words, the OM 
data indicate that, besides the expected synchrotron peak 
likely falling in the infrared band, the SED also presents 
a UV-soft-X-ray bump. 

The existence of a UV-soft-X-ray bump further justi- 
fies the search for a curvature in the X-ray spectrum, and 
the preference given to the double power law model. 
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Table 8. Optical Monitor standard magnitudes of AO 0235+16; numbers in brackets are la errors. 



Date 


V 


B 


U 


UVW^l 


UVM2 


Feb. 10, 2002 
Jan. 18-19, 2004" 
Aug. 02, 2004 
Jan. 28, 2005 


19.22 (0.07) 


20.09 (0.07) 


20.12 (0.10) 
19.60 (0.08) 
20.24 (0.11) 


17.78 (0.06) 
19.95 (0.11) 
19.97 (0.13) 
20.24 (0.13) 


19.68 (0.22) 
19.80 (0.27) 
19.90 (0.20) 



" from Raiteri et al. (i2005i 
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Fig. 11. Ultraviolet to X-ray spectral energy distributions 
of AO 0235+16. The OM optical-UV data have been ob- 
tained by differential photometry with a ~ 1 arcsec aper- 
ture radius to avoid the contamination by ELISA; the pink 
dot refers to February 10, 2002, red triangles to January 
18-19, 2004, green diamonds to August 2, 2004, and blue 
squares to January 28, 2005. The X-ray spectra are shown 
according to Model 5 (double power law), the spectrum 
thickness representing the standard deviation of the data 
from the model fit, and the colours being the same of 
the corresponding OM data. The orange spectrum is from 
Chandra. 

5. Interpretation of the UV-soft- X-ray bump 

The existence of an SED extra-component in the UV- 
soft-X-ray energy range has already been noticed by 
Junkkarinen et al. (200 1 : when analysing the February 
1998 ASCA spectrum together with a contemporaneous 
optical-UV HST spectrum, and by Raiteri et al. (|2005|l . 
when considering the January 2004 XMM-Newton data 
in conjunction with simultaneous optical data acquired 
at the 2.56 m Nordic Optical Telescope in the ambit of 
the WEBT campaign. Even if Junkkarinen et al. I^ESI 
warned that it could arise from a different behaviour of 
extinction in the DLA intervening system with respect to 
the predicted one, we think that the XMM-Newton UV 
and X-ray observations indicate that the existence of a 
real bump must be considered. 

This bump recalls the "big blue bump" found in the 
SEDs of many AGNs, and which is usually ascribed to 
thermal emission from the accretion disc feeding the su- 
permassive black hole, even if accretion disc models gener- 
ally predict a lower flux in the bump tails (see e.g. Krolik 



119991 and Blaes et al. 120011 for the specific case of the 
blazar 3C 273). It is thus interesting to see whether a 
disc model could explain the bump we see in the SED of 
AO 0235+16. By looking at Fig. [nl one can infer that 
this bump would likely peak at log(t/Fi,/ergcm^^ s^^) > 
— 11.5, which implies a peak luminosity vL,^ = And^vFi^ > 
lO^'^ergs"^, if we assume a flat cosmology with Hq — 
Tlkms^^Mpc^^ and Q = 0.27, leading to a source dis- 
tance d ~ 6 X 10^ pc. The peak frequency would likely be in 
the range log(fc'peak/Hz) — 15.5-16, which in the source rest 
frame translates into t'pgak ~ (0.6-2) x 10"'^^ Hz. According 
to the non-LTE disc models computed by Hubeny et al. 
lf2000() the above requirements can be met only in case of 
a rather massive central black hole (^ 10^ M©) and a high 
mass accretion rate (^ 2M0yr~^). 

The lack of optical-UV spectral information during 
the XMM-Newton pointing of February 2002 does not al- 
low to state that the UV-soft- X-ray bump is a strongly 
variable component. In any case, a variable bump could 
be obtained in the accretion disc scenario with changes 
of the model parameters, for instance the accretion rate. 
A more detailed comparison between the observed SEDs 
and the accretion disc models is beyond the aim of the 
present paper; for the moment we can say that we cannot 
rule out that the UV-soft-X-ray bump we observe in the 
AO 0235+16 SEDs is due to thermal emission from an ac- 
cretion disc. The Fe Ka emission line tentatively detected 
in the X-ray spectra would be a signature of this accretion 
disc. 

But other scenarios are also possible. In a recent paper 
Ostorero et al. 1)2004(1 showed how a bump in the UV-soft- 
X-ray region can be obtained in the ambit of a rotating 
helical jet model (Villata & Raiteri [T999|l . by admitting 
a synchrotron contribution to the soft X-ray radiation. 
This contribution would come from the first part of the 
inner emitting jet, while the major synchrotron emission 
at the lower frequencies comes from just outside. This is a 
possible picture for the explanation of the UV-soft-X-ray 
bump suggested by the XMM-Newton observations. We 
could add that, since the "observed" bump seems to be 
more pronounced than that predicted by Ostorero et al. 
|2004), in the framework of the helical jet model the two 
synchrotron components in the SED would come from re- 
gions separated by a sudden variation of opacity. A similar 
interpretation was suggested by Raiteri & Villata l|2()()3|l 
to explain the SED discontinuity of Mkn 501. 

Furthermore, the existence of an additional SED com- 
ponent besides the classical synchrotron and inverse- 
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Compton ones has recently been postulated for BL 
Lacertae. Indeed, when constructing the SED of BL Lac 
during the BeppoSAX observation of October-November 
2000 with radio data and optical data taken by the WEBT 
collaboration (Villata et al. 2002), Ravasio et al. (2003) 
found a misalignment of the soft X-ray spectrum with 
respect to the extrapolation of the optical one. Even if 
according to Bottcher & Reimer 1)2004(1 this misalignment 
may be due to the time averaging performed, it is inter- 
esting here to briefly recall the four different hypotheses 
discussed by Ravasio et al. H2003fl to explain this feature: 

1) a higher than Galactic dust-to-gas ratio towards the 
source (in the BL Lac case there is no intervening system) ; 

2) the detection of bulk Compton emission, i.e. radiation 
produced by inverse-Compton scattering of UV photons 
from the disc and the broad-line region off "cold" elec- 
trons; 3) the presence of two synchrotron emitting regions 
at different distances from the nucleus; 4) the detection of 
a Klein-Nishina effect on the synchrotron spectrum. 

A problem in modelling the absorption towards the 
source, similarly to explanation number 1) for BL Lac, 
cannot be strictly ruled out also for AO 0235-1-16 (see 
Junkkarinen et al. 201) 1 1. Explanation 2), which re- 
quires rather ad-hoc conditions, seems unlikely for AO 
0235-1-16, for which the very high bulk Lorentz factor 
would push the inverse-Compton emission at too high en- 
ergies. Explanation 3) is similar to the one we favour, while 
explanation 4), which is further discussed by Moderski et 
al. H2005|l . is also a possible alternative interpretation. 

6. X-ray light curves 

As mentioned in the Introduction, AO 0235-1-16 has shown 
IDV at various frequencies. During each of the satellite 
pointings presented in this paper, we searched for possible 
rapid variations of the X-ray flux. 

Light curves of the source were extracted from the 
same regions where source spectra were taken, and pre- 
liminarily binned in 1 s intervals. Only bins belonging 
to the good time intervals selected during the filtering of 
high background periods were maintained and a second 
binning of 300 s was performed in order to increase the 
signal to noise ratio. These binned light curves are shown 
in Figs.mirni For the XMM-Newton observations black 
squares, red triangles, and green diamonds refer to the 
MOSl, M0S2, and pn data, respectively. Dotted lines in- 
dicate the average values, with the same colours used for 
the corresponding light curves. 

In order to investigate the existence of a long-term 
trend and, at the same time, the presence of fast varia- 
tions, we verify the goodness of a linear fit. In Table IHl 
we report, for each observation and each instrument, the 
mean count rate < C > of the binned light curve, its stan- 
dard deviation cr, the parameters of a linear fit y — a + b x, 
the corresponding goodness-of-fit, and the probabil- 
ity that the computed fit would have a value of or 
greater. If this probability is greater than 0.1, the fit is 
believable. In the last column we also show the mean frac- 
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5,0x103 1,0x104 1,5x104 a-OxlO" 2,5x10* :3,OxlO-» 
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Fig. 12. X-ray light curve extracted from the Chandra 
observation of August 20-21, 2000; data have been binned 
over 300 s; the horizontal dotted line indicates the average 
value of the dataset 



tional variation, which is a common parameter to char- 
acterize variability fPeterson 12001(1 : /var = ^ cP- — j < 
C >, where 5'^ is the mean square uncertainty of the count 
rates. 

Figure ^| shows the light curves derived from the 
XMM-Newton observation of February 10, 2002. As one 
can see, at least three episodes of remarkable variability 
can be recognized in the light curves at the beginning (only 
MOS detectors), at about half exposure (especially in the 
pn detector, at about ~ 11000 s), and at the end of the ex- 
posure. They imply variations of about 45% in the count 
rates on time scales t =< C > /{\AC/At\) ~ 40 min. The 
variability is confirmed by the high values of /'^ ^-nd null 
probability reported in Table IHl Notice also the decreas- 
ing trend highlighted by the linear fit calculated on the pn 
data (solid green line in the figure). These results imply 
that there are both long-term and fast count rate changes. 
The pn light curve of this XMM-Newton observation has 
been analyzed within several energy bands between 0.3 
keV and 12 keV by Kadler 1(2005(1 . No time lags between 
the soft and hard energy bands could be detected down 
to a limit of ^ 200 s. However, a small but significant 
spectral steepening with increasing flux is found from the 
extraction of separate spectra for the low and high states 
of the source emission. 

On the contrary, no significant variability has been 
found when analysing the X-ray light curves of the other 
pointings, which found the source in a much fainter state. 
This is already visible in Figs. El and 114m 61 since it is 
not possible to recognize variations which are detected by 
all instruments. A confirmation comes from the statisti- 
cal analysis whose results are shown in Table the best- 
fit line is essentially a constant coincident with the aver- 
age value, indicating that there are no long-term trends. 
Moreover, both the low x^/^ values and the probability 
much greater than 0.1 indicate that this constant value is 
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Table 9. Statistics on the Chandra and XMM-Newton hght curves binned every 300 s and fitted with a straight hne 
y — a + b X 



Detector 


< C > 


a a b 




Prob. 


fva.r 




(cts/s) 


(cts/s) 








Chandra: August 20-21, 2000 


ACIS-S 


0.14 


0.02 0.14 -5.36 X 10"** 


1.00 


0.49 


0.00 


XMM-Newton: February 10, 2002 


MOSl 


1.08 


0.15 1.20 -1.40 X 10"'' 


4.26 


0.00 


0.13 


MOS2 


1.11 


0.15 1.23 -1.40 X lO"'^ 


4.11 


0.00 


0.13 


pn 


3.36 


0.38 3.62 -2.81 x 10"^ 


10.52 


0.00 


0.11 


XMM-Newton: January 18-19, 2004 


MOSl 


0.19 


0.02 0.18 6.82 X 10"' 


0.82 


0.89 


0.00 


MOS2 


0.20 


0.03 0.20 1.03 X 10"* 


1.01 


0.46 


0.00 


pn 


0.47 


0.05 0.46 5.90 x lO"'^ 


1.15 


0.16 


0.03 


XMM-Newton: August 2, 2004 


MOSl 


0.19 


0.03 0.18 1.12 X 10"" 


1.01 


0.46 


0.04 


MOS2 


0.19 


0.03 0.19 4.22 X lO"'^ 


0.98 


0.51 


0.00 


pn 


0.49 


0.04 0.50 -1.31 X lO"** 


0.89 


0.66 


0.00 


XMM-Newton: January 28, 2005 


MOSl 


0.09 


0.02 0.08 -1.08 X 10" ' 


1.05 


0.38 


0.04 


MOS2 


0.09 


0.02 0.09 -4.88 x 10"* 


1.23 


0.13 


0.07 


pn 


0.22 


0.02 0.21 1.43 X lO"'' 


0.67 


0.97 


0.00 



XMM-Kcwlon: February 10, 2002 



XMM-Newton: January 18-19, 2004 





Fig. 13. X-ray hght curves extracted from the XMM- 
Newton observation of February 10, 2002; data from pn 
(green diamonds), MOSl (black squares), and M0S2 (red 
triangles) have been binned over 300 s; horizontal dotted 
lines indicate the average values of the datasets; vertical 
lines highlight the main peaks 




5.0xl0'l 



1.5x10-1 
Time [s] 



Fig. 14. X-ray light curves extracted from the XMM- 
Newton observation of January 18-19, 2004; data from pn 
(green diamonds), MOSl (black squares), and M0S2 (red 
triangles) have been binned over 300 s; horizontal dotted 
lines indicate the average values of the datasets; vertical 
lines highlight the main peaks 



a good description of the data. The picture is further con- 
firmed by the mean fractional variation: only in the case of 
the February 2002 observations /var is greater than 10%. 

7. Conclusions 

We have analysed the data obtained during four XMM- 
Newton and one Chandra pointings of the blazar AO 
0235-1-16 in the years 2000-2005. This analysis has care- 



fully taken into account the contribution to both emission 
and absorption from the complex environment at redshift 
z = 0.524. 

We have fitted the X-ray spectra with different mod- 
els (single power law, logarithmic parabola, broken power 
law, curved model, double power law) and have found that 
those implying curvature seem to give a better description 
of the source spectra. The double power law implies pro- 
nounced spectral curvature at soft X-rays. Such curvature 
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XMM-Newton: August 2, 2004 





2.0x103 4.0x10^ 6.0x103 8.0x10= 1.0x104 1.2x10* 

44me [s] 

Fig. 15. X-ray light curves extracted from the XMM- 
Newton observation of August 2, 2004; data from pn 
(green diamonds), MOSl (black squares), and M0S2 (red 
triangles) have been binned over 300 s; horizontal dotted 
lines indicate the average values of the datasets; vertical 
lines highlight the main peaks 

XMM-Newton: January 28. 2005 




Fig. 16. X-ray light curves extracted from the XMM- 
Newton observation of January 28, 2005; data from pn 
(green diamonds), MOSl (black squares), and M0S2 (red 
triangles) have been binned over 300 s; horizontal dotted 
lines indicate the average values of the datasets; vertical 
lines highlight the main peaks 



is not implausible, and favoured in this paper, with re- 
gard to the expected overlapping of the synchrotron (or 
thermal) and inverse-Compton branches of the SED in the 
energy range discussed here (0.3-10 keV). 

We have speculated about the possible existence of a 
Fe Ka emission feature at ^ 2.8 keV in both the Chandra 
and August 2004 XMM-Newton spectra. If this tentative 
detection is real, it may suggest strong gravitational red- 
shift effects from the inner accretion disc. 



The optical~UV data from the Optical Monitor on- 
board XMM-Newton allow one to reconstruct the optical 
to X-ray SED with simultaneous data. The optical-UV 
brightness level is higher when the X-ray one is higher. 
The SED reveals a UV-soft-X-ray bump, which was pre- 
viously noticed by Junkkarinen et al. 1)2004(1 and Raiteri et 
al. (|2005|). If this bump is intrinsic to the source, one pos- 
sible mechanism producing this extra-component may be 
synchrotron emission from an inner region of the jet with 
respect to the region where the lower-energy synchrotron 
component comes from. However, thermal emission from 
an accretion disc, similar to the "big blue bump" observed 
in many AGNs, cannot be ruled out. Investigating whether 
this extra-component is a common feature of the blazar 
class and discovering the mechanism which is responsi- 
ble for it are two intriguing issues which certainly deserve 
further multiwavelength observing efforts. 

We have detected remarkable short-term variability in 
the X-ray light curves derived from the XMM-Newton ob- 
servation of February 10, 2002, when the source was in a 
bright state. We noticed three episodes where the count 
rate changed by about 45% on a time scale of about 40 
min. On the contrary, no significant variability was evi- 
dent in the X-ray light curves corresponding to the other 
observing epochs analysed in this paper. 
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